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In the standard model (SM), the mass and fla- 
vor eigenstates of the B° meson differ. This gives 
rise to particle-antiparticle oscillations which pro- 
ceed in the SM through weak interaction processes, 
and whose phenomenology depends on the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix. The 
time (t) evolution of P° mesons is governed by the 
Schrodinger equation 




with mass matrix M and decay matrix T. The mass dif- 
ference, Am = rriH — tol ~ 2|Mi2|, between the heavy 
and light mass eigenstates, B° H and B° sL , determines the 
frequency of P? oscillations, a quantity precisely mea- 
sured in Ref. [2j. The mean lifetime, r = 2/(Tl + r#), 
is expected to be equal to the mean B° lifetime within 
1% 0| . The decay-width difference, Ar = T L - T H , is 
predicted in the SM to be 0.096 ± 0.039 ps" 1 [3 and 
was first measured by the CDF Collaobarion [5| and, 
recently, by the DO Collaboration with higher precision 
@. It depends on the CP- violating weak phase differ- 
ence between the B®-B® mixing amplitude and the am- 
plitudes of the subsequent P° and B® decays to com- 
mon final states, <p s — arg(— M12/T12), via the rela- 
tion Ar = 2|r 12 |cos(0 s ). While the SM expectation, 
^sm — 4 x io -3 [4]^ is small, contributions from new 
physics processes to B® mixing can lead to a significantly 
different value of the phase, <p s = 0f M + 0f p . The same 
new physics contribution, <fi^ F , would be present in the 
relative phase between mixing and b — > ccs quark transi- 
tions, 2f3 s = 2/3f M - 4>f p , in which the SM contribution 
is defined in terms of CKM matrix elements by = 
&Yg{-V ts V* b /V cs V* b ) » 0.02 4]. Since both 0f M and /?f M 
are significantly smaller than the current experimental 
resolution, we can approximate 2(3 S — — 0^ p = — S . 
Thus the measurement of a sizable value of 2[3 S incon- 
sistent with zero would indicate new physics. In case 
of a non-zero |Ti2 1 , an analysis of time-dependent decay 
rates of B® mesons to two vector mesons becomes sensi- 
tive to the weak phase 2/3 s , even without information on 
the B® flavor at production, because of the interference 
between CP eigenstates. 

In this Letter, we present the measurement of the B® 
meson mean lifetime t and decay-width difference Ar 
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using P° — > J/00 decays followed by J/tp — > and 
— > K + K~ decays. Charge-conjugate modes are im- 
plied throughout this paper. We also extract information 
about the weak phase 2/3 s . The final state is a mixture of 
CP-even and CP-odd states that are distinguished using 
the angular distributions of the decay products. Since 
the P° is a pseudoscalar and the J/ip and are vector 
mesons, the orbital angular momentum between the two 
decay products can have the magnitudes t = 0, 1, or 2. 
The final state is CP-even in S- and D-wave decays and 
CP-odd in P-wave decays. The angular distributions are 
expressed in terms of three angles, 9t, 0t, and ifa, de- 
fined in the transversity basis jjjj]. The angles 0t and 
0t are the polar and azimuthal angles of the in the 
rest frame of the J/0, where the x-axis is defined by the 
momentum direction of the B® and the xy-plane by the 
— > K + K~ decay plane with a positive y component 
of the K + momentum. The angle ipT is the polar angle 
of the K + with respect to the opposite of the P° flight 
direction in the rest frame. 

The data were collected by the CDF II detector at 
the Fermilab Tevatron pp collider between February 2002 
and January 2007, and correspond to an integrated lu- 
minosity of 1.7 fb" 1 . The CDF II detector [1] con- 
sists of a magnetic spectrometer surrounded by electro- 
magnetic and hadronic calorimeters and muon detectors. 
The tracking system is composed of a silicon micro-strip 
detector |9j surrounded by an open-cell drift chamber 
(COT) [lflj. We detect muons in planes of multiwire 
drift chambers and scintillators (llj in the pseudorapid- 
ity range < 1.0. Charged particle identification is 
provided by the time-of-flight system [l2j , complemented 
by the ionization-energy-loss measurement in the COT 
(dE/dx). Events with J/0 — > decays used in this 

analysis were recorded using a dimuon trigger, which re- 
quired two oppositely-charged COT tracks matched to 
muon chamber track segments with a dimuon mass be- 
tween 2.7 and 4.0 GeV/c 2 . 

In the offline analysis, P° — > J/0 decays are recon- 
structed following the procedure described in Ref. [f|. 
We train an artificial neural network (ANN) to separate 
P° decays from the combinatorial background, which is 
the dominant one. We model the signal with simulated 
events and use data from P° mass sidebands (see Fig. [T]) 
to model the combinatorial background. The input vari- 
ables to the ANN are kinematic quantities, vertex fit 
quality parameters, and particle-identification informa- 
tion obtained from the muon system, the time-of-flight 
detector, and the dE/dx measurements. The require- 
ment on the ANN output is selected by maximizing the 
significance S/ \/S + B on data where S (B) is the num- 
ber of signal (background) events in a ±20 MeV/c 2 win- 
dow around the P° mass peak position. The selected 
sample contains about 2500 B® — ► J/0 decays. The 
resulting mass distribution is shown in Fig. [TJ 

To extract r and Ar, we perform an unbinned max- 
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FIG. 1: Invariant J/0 c/> mass distribution with fit projection 
overlaid. The arrows indicate the sideband regions. 



imum likelihood fit with probability density functions 
(PDFs) depending on mass, lifetime, and transversity an- 
gles. 

For the probability density functions (PDFs) of the 
background, we use empirical models with floating fit pa- 
rameters determined from the data. The background has 
a prompt component and a non-prompt component. The 
mass PDF is parametrized by a straight-line function for 
each component. The lifetime distribution is described 
by a delta function at t = for the prompt component, a 
positive exponential for the long-lived non-prompt com- 
ponent, and a negative and positive exponential for mis- 
measured candidates. All lifetime components are con- 
volved with a Gaussian to account for the lifetime reso- 
lution estimated on a candidate-by-candidate basis. Be- 
cause correlations among the three angles are negligible, 
we factorize the angular PDF as a product of polynomials 
in cos 2 (tV), cos(2c/>t) and cos(V't)- The angular distri- 
butions of prompt and non-prompt background events 
agree within uncertainties, and are chosen to be identical 
in the likelihood function. 

For the signal, the mass distribution is described by 
the sum of two Gaussians. The lifetime and the angles 
p = (cos(6*t), 4>t, cos(i/>t)) are correlated for B° signal 
events. The lifetime-angular distribution without accep- 
tance effects is given by 



d 4 P(R,t) 
dpdt 



<x \A \ 2 Mp)T + + \A ll \ 2 f 2 (p)T + 

+ \A ± \ 2 f 3 (p)T_ + |A ||A|||/ 5 (p)cos(<5||)T + 
+ | J 4|||| J 4 ± |/ 4 (p)cos(<5 x -«5||) 

sin(2/? s )(e- rHt -e- rLt )/2 
+ \A \\A ± \f 6 (p) cos(S±) 

sm(2f3 s )(e~ F » t -e- r - t )/2, 



where 

r± 

Hp) 
Hp) 
Hp) 
Hp) 
Hp) 
Hp) 



2 

sin 



cos 

2 



[(1 ± cos(2/3 s ))e^ rit + (1 T cos(2/3 s )) e - r " 4 ] /2, 

"(0 T )(l-sin 2 (0 T )cos 2 ( c /»T)), 
(0 T )(l-sin 2 (0 T )sin 2 ( c /» T )), 
sm*{ip T ) sin 2 (6* T ), 
— sin 2 (^r) sin(2#T) sui(</>t), 
sm(2i/"r) sin 2 (6» T ) sm(24> T )/V2, 
sin(20 T ) sin(26» T ) cos(t/> T )/V2. 

The quantities Aq, A±_ and A|| are the linear polarization 
amplitudes at t = 0, and S± and 5\\ are the strong phases 
of A± and Au relative to Aq, respectively. 

The lifetime-angle distribution is invariant under each 
of the two transformations (2/3 s — > — 2/3 s , S± — > S± + it) 
and (Ar -> -Ar, 2f3 s -> 2j3 s + tt). Because of this 
four-fold ambiguity, this measurement is insensitive to 
the sign of both 2j3 s and Ar. 

The signal lifetime terms are convolved with the same 
Gaussian resolution function used for the background, 
which employs the candidate-by-candidatc lifetime un- 
certainty. To account for different distributions of the 
lifetime uncertainty between signal and background, 
their PDFs are included in the likelihood. These PDFs 
are derived from sideband-subtracted signal events and 
from sideband events, respectively. 

The angular distribution of B® decays described in 
Eq. ([T]) is modified by detector acceptance as well as trig- 
ger and selection efficiencies. This effect is taken into 
account with an acceptance function, e(p), derived from 
simulated B° B — ► J/ip (f> decays. The factor e(p) is de- 
scribed by a three-dimensional histogram with 20 bins in 
each of the angles. 

We consider possible systematic uncertainties due to 
the signal mass model, the lifetime resolution model, the 
0(3%) contamination by B° — ► J/tpK* decays misre- 
constructed and selected as B® candidates not included 
in the background model, the acceptance description, the 
silicon detector alignment, and the model for the angular 
distribution of the background. The largest systematic 
uncertainty for AT is caused by B° mesons reconstructed 
as B® mesons. The largest contributions to the system- 
atic uncertainty on r are the lifetime resolution model 
and the silicon detector alignment. The dominant source 
of systematic uncertainties on the amplitudes is the an- 
gular background model. 

Under the assumption of CP conservation (2f3 s = 0), 
we obtain 



(1) 



T 

Ar 

ill 



A,,' 2 



1.52 ±0.04 ±0.02 ps, 
0.0761S ± 0.006 ps 
0.531 ± 0.020 ±0.007, 
0.239 ± 0.029 ±0.011, 
0.230 ±0.026 ±0.009. 
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FIG. 2: Lifetime distribution with fit projection overlaid. 
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FIG. 3: Angular distributions with fit projection overlaid. 



The first uncertainties are statistical, the second ones sys- 
tematic. We do not quote a point-estimate of the strong 
phase 6\\ because its likelihood profile is non-parabolic 
due to a symmetry point at S\\ = it which makes an uncer- 
tainty estimate unreliable. This analysis is insensitive to 
the second strong phase S± if 2f3 s — (see Eq.([T])). The 
invariant mass, proper decay time, and angular distribu- 
tions with fit projections overlaid are shown in Figs. [T] to 
[3l The measured mean lifetime is compatible with the 
B° lifetime [l3[ as predicted by theory |3j. The mea- 
sured amplitudes are consistent with the ones observed 



FIG. 4: Regions at the 90 % and 95 % confidence level in 
the 2/3 s -AF plane compared with the SM prediction and 
the region allowed in new physics models given by Ar = 
2|ri 2 |cos(0 s ) g. 



in B° — > J/ipK* decays [l4[ as expected under the as- 
sumption of SU(3) flavor symmetry. 

For the constraints on the CP-violating phase, we con- 
struct a 90 (95) % confidence level region in the 2f3 s -AT 
plane using the likelihood-ratio ordering of Feldman and 
Cousins [15j. We choose this method instead of a point 
estimate, because it is not affected by the bias we observe 
in simulated experiments. The bias is of the order of the 
statistical uncertainty for input values of Ar or 2/3 s close 
to zero, which are near to the SM expectation. The bias 
can be understood from Eq. ([I]). If 2/3 s approaches zero, 
the two terms proportional to sin(2/3 s ) vanish, and this 
analysis becomes insensitive to 8±. The same effective 
loss of degrees of freedom in the fit occurs when Ar ap- 
proaches zero and multiple degenerate solutions for 2/3 s 
and S± exist. 

To obtain the likelihood ratio distribution for given val- 
ues of Ar and 2/3 s , we use experiments simulated with 
values for all other parameters determined by a fit to data 
[l^ l . We checked that alternate choices of these values do 
not affect the coverage properties of our algorithm. Sys- 
tematic uncertainties are not included in the algorithm, 
since they are all negligible. 

The resulting confidence region is shown in Fig. 3] 
Since both B® mass eigenstates have the same angu- 
lar distribution at 2(3 s = ±7r/2, the sensitivity on Ar 
decreases towards this value. For the SM expectation 
(Ar w 0.1 ps" 1 and 2f3 s w 0), we find the probability 
to get an equal or greater likelihood ratio than the one 
observed in data to be p — 22%, corresponding to an 
agreement at 1.2 Gaussian standard deviations. 

In summary, we report the measurement of the mean 
lifetime, the width difference, and the amplitudes in 
B® — > J/0 <j> decays assuming CP conservation. This 
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measurement improves the precision of the current best 
measurement Q by 30-50%. It is in good agreement 
with previous results and the SM expectation. In addi- 
tion we derive constraints on Ar and the CP-violating 
phase 2(3 S . Our data are consistent with the SM expec- 
tation of 2f3 s » 0, but sizeable values allowed within new 
physics models cannot be ruled out. 
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